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A covariant Fokker-Planck type equation for a simple gas and an equation for the Brownian 
motion are derived from a relativistic kinetic theory based on the Boltzmann equation. For the 
simple gas the dynamic friction four-vector and the diffusion tensor are identified and written in 
terms of integrals which take into account the collision processes. In the case of Brownian motion, 
the Brownian particles are considered as non-relativistic whereas the background gas behaves as a 
relativistic gas. A general expression for the semi-relativistic viscous friction coefficient is obtained 
and the particular case of constant differential cross-section is analyzed for which the non-relativistic 
and ultra relativistic limiting cases are calculated. 
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I. INTRODUCTION 

The problem of Brownian motion played a fundamen- 
tal role in the early verifications of kinetie theory after 
the estimation of the Avogadro number by Perrin in 1908 
on the basis of Einstein's theory. Brownian motion and 
diffusion processes have a wide range of applications that 
goes from chemistry, solid state, quantum physics, plas- 
mas, astronomy and astrophysics to social sciences, life 
sciences and biology 

Other kinds of applications can be found if we extend 
the Brownian motion to the regime of the theory of rel- 
ativity. The relativistic Brownian motion could have im- 
portant applications in plasma physics, in high energy 
physics [31, in astrophysics - for example in the analysis 
of the gamma ray burst jets where a Brownian motion 
of the electrons in the electrostatic wave field could be 
present Q ~ in relativistic corrections to the Sunyaev- 
Zeldovieh effect and also in the evolution of dark 
matter Q. 

Recently, several different approaches of a relativistic 
theory of stochastic processes have been used for the im- 
plementation of the Brownian motion concept into the 
theory of special relativity @, 0, [l] , some of them gener- 
alized their results to a curved space-time , and studied 
the problem from the point of view of reaching quantum 
theory with Nelson's method [l^l , and also from quantum 
gravity [ill , [l^ . A comparison among different relativis- 
tic stochastic processes existing in the literature can be 
found in p^ . 

In the works [ij, [l^ , Dunkel and Hanggi constructed a 
relativistic theory for the Brownian motion starting with 
a relativistic version of the Langevin equation and ob- 
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tained that the related Fokker-Planck equation can be 
written as a continuity equation, but with different ex- 
pressions for the flux density by using different interpre- 
tations of the stochastic processes. Moreover, they found 
that only one of the approaches leads to the equilibrium 
relativistic MaxwcU-Jiittner distribution function, while 
the others differ through factors that depend on the en- 
ergy of the particles. In a recent paper by the same au- 
thors [iBl they have started from a microscopic collision 
model and constructed a criteria to identify the equilib- 
rium distribution of the particles. Firstly they have used 
this criteria in the non-relativistic case and recognized 
the Maxwellian distribution function, then they applied 
the same criteria in the relativistic case and, following 
their previous results, found that the relativistic distri- 
bution function that satisfies their criteria, differs from 
the relativistic Maxwell-Jiittner distribution function by 
a factor proportional to the inverse of the relativistic ki- 
netic energy. This result does not agree with the one 
which comes out from the relativistic kinetic theory based 
on the Boltzmann equation where the only distribution 
function which implies a vanishing collision term in equi- 
librium is the Maxwell-Jiittner distribution function (see 
e.g. [13, HI and [3). 

It is important to have a theory of Brownian motion 
founded on kinetie theory, since that theory could give 
the connection between the fundamental microscopic dy- 
namics and macroscopic schemes which have measur- 
able properties. With such a theory one can obtain the 
equations that govern the motion, and also calculate the 
transport coefficients as a function of the properties of the 
specific system. Is the aim of the present work to inves- 
tigate a relativistic generalization of the Fokker-Planck 
equation and of the equation for the Brownian motion 
within the framework of the relativistic kinetic theory 
based on the Boltzmann equation. 

This work is structured as follows. In Sec. II a co- 
variant version of the Fokker-Planck equation is derived 
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from a relativistic kinetic theory based on the Boltz- 
mann equation and apphed to a system with only one 
constituent. The same assumptions of grazing collisions 
valid for a non-relativistic gas (see e.g. |20j and [2l|) are 
also considered here, which consist in small deflections on 
the scattering angle and small changes in the momentum 
of the particles at collision. The dynamic friction four- 
vector and the diffusion tensor arc identified and written 
in terms of integrals which take into account the collision 
processes. In Sec. Ill a mixture of two species is consid- 
ered where one of the components has a small particle 
number density and whose particles have a large mass 
when compared with the other component. Those con- 
ditions allow us to call this case a Brownian motion by 
analogy with the fluctuation process (see e.g. [Ill)- As in 
the non-rclativistic case (sec e.g. the former is iden- 
tified with the Brownian particles and the latter with the 
background gas. The Brownian particles are considered 
as non-relativistic whereas the background gas behaves 
as a relativistic gas which may alternate from the non- 
relativistic limit to the ultra-relativistic limit depending 
on the ratio between the rest energy of the particles of 
the gas and of the thermal energy of the mixture. A gen- 
eral expression for the semi-relativistic viscous friction 
coefficient is obtained and the particular case of constant 
differential cross-section is analyzed for which the non- 
relativistic and ultra relativistic limiting cases are calcu- 
lated. Concluding remarks are given in Sec. IV. 



II. RELATIVISTIC FOKKER-PLANCK 
EQUATION FOR A SIMPLE GAS 

Let a particle of a simple relativistic gas be character- 
ized by its rest mass m and by the space-time coordinates 
(x") = (c<,x) and momentum four-vector {p") = (p°,p), 
where the component of the momentum four-vector p'^ 
is constrained by p° = ^/\p\+rn^ . The state of the 
relativistic gas in the phase space is characterized by the 
one-particle distribution function /(a;",p") = /(x, p, i), 
such that f{x,p,t)d^x(fip gives at time t the number of 
particles in the volume element (Px about x and with 
momenta in a range cfip about p. 

For an elastic collision of two particles with momen- 
tum four- vectors denoted by p" and p", the energy- 
momentum conservation law reads p" + p'^ = p'" 



where the quantities p'", p'^ denote the values of the 
momentum four-vectors after collision. 

The one-particle distribution function satisfies the re- 
lativistic Boltzmann equation (see e.g (isl. [Toj) 
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df_ 



dfK°' 
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if'J'-fJ) Fadn 



d^p^. 
P*o 



(1) 



where the usual abbreviations = /(x, p'^,t), /' = 
/(x,p',t), /, = /(x,p,,i), and / = /(x,p,<) were in- 
troduced. In the above equation, a denotes a differential 
cross-section, dfl an element of solid angle which charac- 
terizes the scattering process, K" the Minkowski external 



force and F = \J {p"Pa)^ — m'^c^ the so-called invariant 
flux. 

Under the assumption of grazing collisions that could 
take place in long range interactions, only small changes 
in the momentum of the particles occur due to small 
deflections in scattering angle. Then the collision term 
in the right-hand side of the Boltzmann equation p]). 
denoted by Q(/, /*), can be approximated as follows. 

First new variables are introduced, namely, the total 
P" and the relative momentum four-vectors, defined 
by 



P" =p"+p: = P'", Q°'^p°'-p°, 

(2) 

For these quantities the following relationships hold 
P"Qa = and = + Am'^c^ where P^P^, = P^ and 

Further, from the energy-momentum conservation law 
and eq. ^ one can write the differences between the 
post- and pre-coUision momentum four-vectors as 




-^pT 



p. 



i(0'"-Q") = -iAQ". 



(3) 



Hence for small changes of the momentum of the par- 
ticles at collision one can expand the one-particle distri- 
bution function in Taylor series, which up to the second- 
order terms reads 

/(p")-/(p-) + iAQ'§ + iiQ-A(,^^. (4) 

with a similar expression for f(j>'i), by making the 
changes p^ — > and AQ* —AQ^. In these expres- 
sions we take only the first two derivatives because the 
other contributions are of smaller order thanks to the 
hypothesis of grazing collisions. 

Now it is possible to approximate the collision term of 
the Boltzmann equation ([1]) as 



Q(/,/*) = J 



(//* 



Fadn 



d^p^ 
P*o 



with the help of the relationship 



d 



_d_ 

dpi 



_d_ 

dp^ 



(5) 



(6) 



In order to transform the integral (O, the center-of- 
mass system is chosen where the spatial components of 
the total momentum four-vector vanish, i.e. 



(P") = (P",0), 



(Q") = (0,Q). 



(7) 



Now the element of solid angle in ([5|) can be written as 
dft = sinQdQd^, where 8 and <!> are polar angles of 
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Q'" with respect to and such that represents the 
scattering angle. Further, without loss of generality, 
is chosen in the direction of the 3-axis, so that one can 
write Q°' and Q'" as: 






sin O cos $ 
sin O sin $ 
cos 



(8) 



By using the above representations it is easy to calcu- 
late the following integrals in the variable < $ < 27r, 
yielding 



Ag*— cr sine de 

P*o 



(9) 



By invoking the divergence theorem again and using 
it follows that 



(14) 



Q(/,/*) 



d_ 



dJD'- 



dpi 



(15) 



where the spatial components of the coefficient of dy- 
namic friction and the diffusion coefficient Z?'^ are 
given by 



A' = / f^AplFadfl 



d^p^ 
P*o 



(16) 



/ AQ'AQ^— a sine dOd^ = -(QV^ + Q'Q^)Q^- 
J P*o 

(10) 

Note that that the differential cross-section is a function 
of cr = cr((5, Q) whereas the invariant flux is given by 
F/p*o = Q- Above S denotes the following integral 



E = 27r y"(l - cos e)cr sin 9 dQ, 



(11) 



and 77*-' are the spatial components of the metric tensor 
(j]ap) = diag(l, — 1, — 1, — 1). In order to obtain the in- 
tegral ()10p the following approximation was taken into 
account sin^ w 2(1 — cos 9), since only grazing colli- 
sions between the particles with small scattering angles 
are considered. 

By differentiating pU]) with respect to and consider- 
ing the relationship dQ/dQ^ = —Qj/Q, one can obtain 
the following connection between the integrals ([9]) and 



d 



AQ'AQ^Qadn = 2 / AQ'Qadn 



(12) 



Now by making use of equations ^ and , the col- 
lision integral ([5]) becomes 



a(/-/.) = i/|i{ 



m^AQ'AQ^Qadn] d^p. 



l_d_ 



djff. 



■AQ'AQ^Qadfld^p^ 



(13) 



The first term on the right-hand side of the above equa- 
tion vanishes, since the hypothesis of grazing collisions is 
used and it is possible to convert - thanks to the diver- 
gence theorem - the volume integral in the momentum 
space into an integral at an infinitely far surface where 
the distribution functions tend to zero. The second term 
on the right-hand side can be manipulated by using equa- 
tions ([6]) and ((12]), yielding 



2 / fJ^AQ'QadVLd^p^ 



D'^ = i / f^AplApi Fadn^ 
2 J p*o 



(17) 



It is important to call attention to the fact that left- 
hand side of the Boltzmann equation is a scalar invariant, 
but the collision term written in the expression (jl5p is 
not an invariant. In order to write a covariant form of 
the Boltzmann equation one has to recall that the above 
results were obtained by considering the center-of-mass 
system. In this system Ap^ = 0, and one can include the 
zcroth components 



d^p^ 
P*o 



D^o = D°^^U UApMFadil^, 
2 J p*o 



D°'^U UAplAplFadn^ 
2 J p*o 



(18) 



(19) 



(20) 



into the collision term p5l) so that it becomes a scalar 
invariant. While the components of the coefhcient of dy- 
namic friction and of the diffusion coefficient are 
related with the changes of the momenta Apl at collision, 
the components A" and D^^ take into account the energy 
changes Ap" and the changes of both Ap" and Apl . 
They are the temporal components of the four-vector A" 
and of the four-tensor D"^. 

Hence the relativistic Boltzmann equation ([T]) reduces 
to the relativistic Fokker-Planck equation, namely. 



P 



dl_ 



d 
dp^ 



a/3 



dpP 



(21) 



For a spatially homogeneous case without external 
forces, the Fokker-Planck equation (PT|) reduces to 



/a/ 

c dt 



(22) 
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which represents a continuity equation in momentum 
space with the particle flux density J^" given by 



III. RELATIVISTIC EQUATION FOR 
BROWNIAN MOTION 



(23) 



The stationary solution of cq. ([221) for the distribution 
function / is obtained by assuming that the particle flux 
density vanishes, i.e. JF" ~ 0. This condition is equiv- 
alent to the one that in equilibrium the collision term 
of the Boltzmann equation vanishes and the distribution 
function is characterized by the Maxwell- Jiittner distri- 
bution [H, dll, namely, 



/ oc exp 



kT 



(24) 



where Ua - such that U^Ua = - is the hydrodynamical 
four- velocity of the gas. 

If we insert into by considering = 0, it 

follows that 



kTA" = -UfiD"^, where A" = A" 



(25) 



The representations of the four- vector A" and of the four- 
tensor D'^P in terms of the four-momentum read 



(26) 



where A, D and D are scalar coefficients and 13" = 
-4£> - D. 

The insertion of the representations (pS)) into the equa- 
tion p5l) i leads to 



D = 0, kTA^^^D, 



(27) 



i.e., only one among the three coefficients is linearly in- 
dependent. 

Hence, thanks to and ^7\ . the spatially homoge- 
neous Fokker-Planck equation p2|) reduces to 



jf_dl _d_ 
~c'dt ^ dp^ 



Dp'' f p'^Ug 
■' kT 



0. (28) 



In a Lorentz rest frame - where = (c, 0) - one can 
obtain from equation (|27p 9 



ikTA 



D 



V 

2^ 



(29) 



where the leading term, for v <^ c, corresponds to the 
usual relation between the diffusion coefficient and the 
friction coefficient. 



In this section a mixture of two constituents is con- 
sidered, where one of the components consist of heavy 
particles of rest mass mb while the other by light parti- 
cles of rest mass mg, so that mb 3> rug. The component 
with light particles describes a rarefied gas with particle 
number density much larger than that of the constituent 
with heavy particles and which characterizes the Brow- 
nian particles rif, <^ Ug. The gas is supposed to be at 
equilibrium with a Maxwcll-Jiittner distribution function 



_ f(0) _ 



/g — ft 



4TrmlckTK2iCg) 



(30) 



Above, rig denotes the particle number density of the 
gas, c the speed of light, k the Boltzmann constant, 
whereas T and are the temperature and the four- 
velocity of the mixture, respectively. The symbol i^2(Cg) 
refers to a modified Bessel function of second kind and 
Cg = uigC^/kT gives the ratio between the rest and ther- 
mal energies of the gas particles. 

The two assumptions above for the Brownian con- 
stituent - that it has a small particle number density 
and particles with large mass with respect to the gas 
constituent - imply that it may be considered as a non- 
relativistic gas with negligible collision term with respect 
to its particles. Hence the Boltzmann equation for the 
one-particle distribution function of the Brownian parti- 
cles can be written as 



5/b , 9/b 9(/bF') 

dt ^^dx^ dpi 



(/i°^7b-/i"Vb)50'TdM3pg. 

(31) 

The above equation follows easily from the Boltzmann 
equation ([T]) written for a Brownian particle through its 
multiplication by c/p^, the introduction of the velocity 
^b = cpb/Pb ^^'^ of the M0ller velocity defined by 



90 



cF 



(vg 



Vb 



r(vg X vb)2 



(32) 



Since the collisions of the gas particles with the Brown- 
ian particles affect very little the depart from equilibrium 
of the latter, one can suppose that the distribution func- 
tion of the Brownian particles can be written as 



(33) 



where fj^^^ refers to the exponential factor of a Maxwell- 
Jiittner distribution function for the Brownian particles 
and h{pU represents a deviation from this distribution 
when we assume space homogeneity. Moreover, based 
on this assumption one can expand the deviation for the 
post-collision momentum h{p'^) in Taylor series and by 
keeping up to second order terms in the difference Ap^ ~ 
P'b - Pb' yields 



dh 



1 



/^bb) = HpD + Apl— + -A^Ap, 



' dpldpl^ 



(34) 
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Equations ((33|) and p4|) are introduced into the colli- 
sion term of the Boltzmann equation pip so that it can 
be written as 



Q(/b,/g)=/r-4x 



(35) 



where |gp = g'^g^. For the second integration the ele- 
ment (Ppg = sm(j)d(j)d'&\pg\'^d\pg\ is written in terms of 
the spherical coordinates (|pg|, 0), where and (j) are 
the spherical angles of with respect to p]^. Hence the 
integrations of and g^g^ with respect to i9 become 



where the relationship /^^°Vf'^°^ 



Above A 
integral 



~ Jh J 

Ug/ A'Km'ickT K2{Cg) and J is the following 



^/g°^ was used. 



2tt 



g'dd = 27rc 



77° 



cos(/.-l)^, (42) 
iPbl / Pi 



~ g^ad^d Pg 



(36) 

The integral ((36)) can be transformed as follows. First 
one can note that for a non-relativistic Brownian particle 



Ph — Ph so that energy conservation law leads to p2- 



p'°- 



Further one can introduce a relative velocity defined by 



CPg 

Pi 



CPb 

pI'^ 



CPg 

Pi 



pI' 



(37) 



such that the difference Ag* = g'* — (7* can be written, 
thanks to the momentum conservation law, in terms of 
the difference Ap* as follows 



g^g^dd = 2t:c^ 



1 IPgP [pl? .2 

2 [pl? IPbP ' 



P° IPbl 



cos ( 



p\A 



1 Ipg 



(43) 

The invariant flux F — ^ (pgap^)^ — m^rn^c^ can be 

approximated for the case of a non-rclativistic Brownian 
particle by 



pIp190 0, I ( 1 IPb| Pg , 

^.pjp,| 1-^^COS0 



(44) 



c 



(38) 



For a relativistic gas of rest massless particles |Pg/Pg| = 
1 whereas |Pb/Pb| = Wf,/c, hence one can approximate the 
difference given by ([38]) as 



(39) 



« --,AK. 



The above approximation is also valid for a gas with mas- 
sive particles, since mb 3> Wg and Vg ^ c imply that 

The integral (|36p can be written in terms of the dif- 
ference Ag' by using ([39)) . In the resulting equation, 
two integrations can be performed, the first one can be 
done by introducing the scattering angle x and the az- 
imuthal angle e which are the spherical angles of _g'* with 
respect to g*. By writing the element of solid angle as 
dfi = svaxdxde- and integrating the differences Ay* and 
Ag^Ag^ with respect to the azimuthal angle < e < 27r, 
yields 



by considering terms up to the first order of |pb|/Pb = 
fb/c. Note that the above approximation is valid if 
Pg/|Pg| does not diverge, but in the following wc shall 
show that after some transformations such a term disap- 
pears from the final form of the integral X. 

Since the differential cross-section is a function of the 
invariant fiux and of the scattering angle it can be ap- 
proximated by 

a = a{F,x) ~ fT(|PgU) [I - ^^-^p^cos, 

(45) 

Now by using ([55)) through ([i5)) in equation ([M)) and 
integrating the resulting equation with respect to the an- 
gle < (f> < TT it follows that 



|PgP(7(|Pg|,x)(l-cosx)sinx 



p{ 
dpi pI 



l + |Pg 



din a 
5W 



27r 



A5Me = -27r(l-cosx) g\ 



(40) 



dp\M' (Pg)' 



dx 



c^lPg 



(46) 



Ag'Ay^de = 27r| (1 - cosx) 



In order to integrate by parts the derivative of the cross 
section a with respect to |pg| a Lorentz rest frame in 
which (J7") = (c, 0) is chosen, so that one can obtain 



■ sin^ X 



(41) 



; da 
5|Pg 



■d|Pg 



4|PgP 
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kTp' 



(^(|Pg|,x)t^|Pg|- 



(47) 



Now the substitution of the above resuh into leads 
to 

8 f '=pg 
I = ---k'^c j e-"cr |pg|5cr(|pg|,x)(l - cosx)sinx 



dh 



V 



dpldpl 



Pi 



(48) 



The final form of the collision term (|35p in terms of 
the distribution function of the Brownian particle can 
be obtained from ((48|) by using the representation ([33ll . 
where the differentiation was also performed in a Lorcntz 
rest reference frame, namely 



Q(/b,/g) = ^ 



dpldpi dpl 



(49) 



so that the Boltzmann equation for the Brownian particle 
can be written as 



dt dpi 



dp\M dpl 



(50) 



In and ([50)1 77 is the so-called viscous friction coeffi- 
cient for a scmi-relativistic case and it is given by 



7/ 



3 mb{mgkTyK2iCg) 



I CT(|pg|,x)(l-cosx)sinxdxe |Pg|^^^ 

J Pe 



(51) 



Equation ()50|1 is the extension to a background gas 
of relativistic particles of the Fokker-Planck type equa- 
tion for Brownian motion which was analyzed by Chan- 
drasckhar [24|, Green and Wang Chang and Uhlen- 
bcck for the case of a non-relativistic background 
gas. The relativistic corrections appear here in the vis- 
cous friction coefficient. 

In order to evaluate the integral for the viscous friction 
coefficient (|5ip it is necessary to know the differential 
cross-section a. For a constant differential cross-section 
the integration in the angle < x ^ '''' is straightfor- 
ward and by introducing a new variable y — Pg/rUgC so 

that |pg| = mgCi/ — 1, the friction viscous coefficient 
becomes 



3 nihimgkTy K2{Cg) 



-CeV 



{y'-lfdy, (52) 



which leads through integration to the final form of the 
viscous friction coefficient of the Brownian particles in a 
background gas of relativistic particles: 



327rngCrfcT e-«« 
3mhcK2{Cg) Q 



(3 + 3Cg + CD 



(53) 



For low temperatures the rest energy of the background 
gas irigC^ is much larger than the thermal energy kT, 
so that (g ^ 1 and the semi-relativistic viscous friction 
coefficient can be approximated by 



^^J2mgTTkT 
3mb 



1 
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128(1 



(54) 



By using the hard-sphere cross-section of a non- 
relativistic gas a = (P /4 - where d denotes the diam- 
eter of a particle - the first term in the expression (|54|) 
reduces to the result obtained by Wang Chang and Uh- 
lenbeck [2^. The other terms in the series are related 
with relativistic corrections. 

For very high temperatures the parameter (g <^ I and 
this condition characterizes the ultra-relativistic regime. 
In this case the approximation for the viscous friction 
coefficient reads 



IGTTTlgCrfcT 

nihC 




(55) 

where 7 = 0.577215664. . . is the Euler constant. In this 
case the first term is the leading one and the others are 
corrections to the former. 



IV. CONCLUDING REMARKS 

In this work we obtain a manifestly covariant relativis- 
tic Fokker-Planck type equation for the evolution of the 
distribution function of a simple gas under the assump- 
tion of grazing collisions. We also get a semi-relativistic 
Fokker-Planck equation for the Brownian motion. 

Although the Fokker-Planck equation pT|) and the 
equation of the Brownian motion (j50p arc alike - at least 
in the non-relativistic limit - one has to be very careful 
with the temptation to obtain the latter from the for- 
mer. One can write a Fokker-Plack equation like the 
one given in (j2ip for a mixture of two constituents where 
one of constituents has a small particle number density 
with respect to the other, so that only collisions between 
dissimilar particles are taken into account. However, in 
Sec. II only grazing collisions between the particles were 
taken into account, while in the case analyzed in Sec. Ill 
the collisions between the particles of the background gas 
and the Brownian particles were not restricted to grazing 
collisions. 

The extension of the Fokker-Planck type equation for 
a simple gas in the presence of a gravitational field is 
straightforward, since the Boltzmann equation in gravi- 
tational fields is written as (see e.g. [l^) 



^ op" J p*0 

(56) 

Here ^ = ^ ~ det(gpy) where g^^^ denotes the metric 
tensor and FJ^^^ is the Cristoffcl symbol. Indeed, the 
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Fokker-Planck equation for this case becomes 



df_ 

dxf^ 



dp"' 



d 



dfD 



}1U 



dp" 



(57) 



where the dynamic friction four-vector and the diffu- 
sion tensor D^'^ are given by 



J P*Q 



(58) 



Recently in [2^ is shown that a modified Maxwell- 
Jiittner distribution could be obtained from a modified 
principle of maximum entropy that take into account the 
properties of Lorentz group. This is a hint that the rela- 
tivistic Boltzmann equation might require a slight modi- 
fication. Indeed the relativistic kinetic theory considered 
here is not the only version, there are others proposals 
like [131 and [1^ which try to improve some points of the 
theory. 



= i / f^Ap^Ap'CFadn^^. (59) 
^ J P*o 

In Sec. Ill we made an analogy with Brownian motion 
in the relativistic case, that analogy was made through 
scattering assumptions based on specific properties like 
mass and density ratios of the constituents of the sys- 
tem. These assumptions allow us to made a physical 
analogy with Brownian motion and none of the math- 
ematical stochastic properties of Brownian motion were 
used. In fact the evolution in this version of Brownian 
motion is provided by the collisions while in other works 
@i 01 @i [3 and [iBl is given by a stochastic term. 
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